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I. INTRODUCTION 



Direct detection of Dark Matter (DM) relies on signals due to energy deposited from DM 
recoiling on matter in a detector. The DAMA collaboration has provided strong evidence 
for an annually modulated signal in the scintillation light from sodium iodine detectors. 
The combined data from DAMA/Nal [1] (7 annual cycles) and DAMA/LIBRA [2] (4 annual 
cycles) with a total exposure of 0.82 tonyrs shows a modulation signal with 8.2 a significance. 
The phase of this modulation agrees with the assumption that the signal is due to the 
scattering of Weakly Interacting Massive Particles (WIMPs) forming the DM halo of our 
Galaxy. 

After a series of experiments with negative results, the CDMS collaboration has re- 
cently presented an analysis where 2 events over an expected background of 0.9 ± 0.2 events 
have been seen after 612 kg days of exposure [3]. The probability to observe two or more 
background events is 23%, which means that the two events neither provide a statistically 
significant evidence for DM interactions nor can they be rejected as a background. In this 
paper we investigate both hypotheses and confront them with the results of the other DM 
direct detection experiments. Another recent direct detection result is due to XENONIO [4] 
that is searching for signs of DM scattering on liquid xenon. In a recent re-analysis [5] of 
their 316.4 kg day exposure they found 13 events with expected background of 7.4 events, 
which was used to set limits on DM scattering cross sections. 

Both of these very recent experimental results merit re-examination of the DAMA signal 
and its consistency with the results from the other direct detection experiments. In the 
present manuscript we address the following questions: 

• Is it possible to reconcile the DAMA annual modulation signal with the constraints 
from all other experiments? 

• If interpreted as DM signal, are the two events observed in CDMS consistent with 
results of other experiments? 

We will focus on four classes of DM scattering cross sections that cover a large set of DM 
models: the spin-dependent (SD) and spin-independent (SI) WIMP-nucleon scattering that 
can be either elastic or inelastic. We will use the shorthand notations eSD, eSI, iSD, iSI, to 
denote these four classes. 

Part of the answer to the first question is already well known from the literature, since 
many of the interpretations of the DAMA signal are in conflict with the constraints from 
other DM direct detection experiments. For instance, recent analyses of the eSI case have 
been performed in [6-14]. We will show that this scenario gets even more disfavoured due 
to the new CDMS and XENONIO results. Similarly, the explanation of the DAMA signal 
due to the elastic spin-dependent scattering (eSD) is disfavoured by the strong constraints 
from COUPP [15], KIMS [16], and PICASSO [17]. 

The case of inelastic spin-independent (iSI) scattering of a DM particle to a nearly de- 
generate excited state [18] is also tightly constrained, especially by CRESST-II [19] and 
ZEPLIN-III [20], c.f. recent analyses [21-26]. However, as we will show, the spin-dependent 
inelastic scattering (iSD) offers a viable explanation of the DAMA signal, consistent with 
all other constraints. To the best of our knowledge, this solution to DAMA has not been 
discussed in the literature before. 
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Unlike DAMA the CDMS signal is statistically very weak. Even so, focusing on the 
second question above, we also entertain the hypothesis that the two events in CDMS are 
indeed due to DM and perform a maximum likelihood fit allowing for the presence of a 
signal. In the case of elastic scattering (eSI and eSD) we obtained a very weak (at the la 
level) indication for a positive DM signal, and compare the corresponding allowed regions 
to the constraints from all other experiments. 

The outline of the paper is as follows. In section II we collect expressions for predicted 
event rates in DM scattering experiments, in section III we describe details of individual 
experiments relevant for the combined analysis, while in section IV we then give the re- 
sults of the fits. A simple DM toy model leading to inelastic spin-dependent scattering is 
presented in section V, followed by conclusions in section VI. In appendix we also include 
the interpretation of most recent CoGeNT result that appeared after the publication of the 
manuscript. 

In addition to the four general classes of DM models considered in this work, there 
are also other proposals to explain the DAMA signal, including for example mirror world 
DM [27], DM with electric or magnetic dipole moments [28], leptophilic DM [29], resonant 
DM-nucleus scattering [30], modified scattering due to DM form factors [31, 32], and atomic 
DM [33]. The implications of the new CDMS results for the above models is an interesting 
open question that, however, is beyond the scope of this paper. 

II. EVENT RATES 

The differential counting rate in a direct DM detection experiment (in units of counts 
per energy per kg detector mass per day) is given by 



where Ed is the energy deposited in the detector, po is the local DM density (which we take 
to be 0.3 GeVcm"^), r] is the number density of target particles, pdet is the mass density of 
the detector and da/dEd is the differential cross section for scattering on a target nucleus. 
If the target contains different elements (like in the case of the DAMA Nal crystals), the 
sum over the corresponding counting rates is implied. 

In eq. (1), /©(v) is the local WIMP velocity distribution in the rest frame of the detector, 
normalized according to J d^v fQ{\) = 1. It follows from the DM velocity distribution in 
the rest frame of the galaxy, /gai(v), by a Galilean transformation with the velocity of the 
Sun in the galaxy and the motion of the Earth around the Sun. For /gai(v) we assume the 
conventional Maxwellian distribution with v = 220kms~^ and a cut-off due to the escape 
velocity from the galaxy of fesc = 650kms~^: /gai(v) oc exp(— v^/t;^) — exp{v'^^^/v'^) for 
V < Vesc and zero for v > Vesc- The precise value of the escape velocity has a negligible 
impact on our results for elastic scattering, while in the inelastic case the precise shape of 
the tails of the velocity distribution is important [34], see the discussion in sec. IV. 

The lower limit of the integration in eq. (1) is set by the minimal velocity timin that the 
incoming DM particle has to have in order to be able to deposit an energy Ed in the detector. 




(1) 
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For the case of inelastic ~^ x'^ scattering it is given by 

/ niNEd 



+ ^]. (2) 



where /x^^a? = m^rriN / {m^ + itin) is the reduced mass of the nucleus-DM system, with niN 
and the nucleus and DM masses respectively, while S is the mass difference between x' 
and X- The same equation also applies to elastic x^ ~^ scattering, with 5 = 0. As 
observed in ref. [18] for appropriately chosen S one can suppress the signal in experiments 
where DM scatters on lighter nuclei, while not significantly affecting the rate in DAMA 
(see also [35, 36]). Namely for 6 3> mjqEd/ ij^xN the minimal velocity v^^^ falls with mjq. 
If the signal is coming from the tails of the velocity distributions, the difference between 
lighter and heavier nuclei, such as germanium vs. iodine, can be significant (for v^^i^ > v^sc 
the scattering is completely absent). Furthermore, the inelasticity also suppresses the low 
energy signal, changing the shape of the expected event rate from an exponentially falling 
function of the recoil energy to a bump-like signal at higher energies. This, in addition, 
improves the fit to the DAMA modulated signal energy spectrum. 

The differential cross section for scattering on a target nucleus is (per assumption) given 
by the spin independent (SI) and spin dependent (SD) contributions, which are convention- 
ally written as (see e.g. [37]) 

where o"^^'^^ are the integrated SI and SD cross sections for DM scattering on nucleus, 
but with form factors factored out. For the SI form factor F{Ed) we use [38] F{Ed) = 
3g-«; s /2|^gj]^('^j.^ _ cos(/€r)]/(Kr)^, with s = 1 fm, r = y/K^ — 5s^, R = 1.2A^^^ fm, 
K = y/2rnNEd (and ~ — k^). The SD form factor S{Ed) is computed according to ref. [39] 
for ^'^^Cs (abundant in the Csl crystals used by the KIMS experiment) and according to 
ref. [40] for all other nuclei. 

Even though the form factors were factored out of the definitions of o"^^'^^, these quan- 
tities still depend on nuclear structure through isospin content (the number of protons vs. 
neutrons). The SI cross section is thus 



SI ^ l^fp + ~ ^)/n]^ ^^XN SI (A\ 
^ f2 ,.2 ' 



with A the atomic mass number, Z the charge of the nucleus, the SI DM couplings to 
proton and neutron respectivelly, fi^p the reduced DM-proton mass, and (7^^ the SI cross 
section for scattering of DM on a proton. In the fits we will assume fp = fn for definiteness 
and quote results in terms of 0"^^. Since the ratio A/Z is similar for different nuclei this 
choice mostly affects only the overall value of a^^ , while it does not affect the relative sizes 
of contributions from different experiments. It is easy to rescale our results for different 
values of and /„ through af ^ cr^'/{Z/A + (1 - Z/A)fjfpf. 

The SD cross section depends in addition on the spin J of the nucleus 

^'°^(^'^) = o 2^2^o^LiJ ^o^oo(g) + aoai5oi(g) + a?5n(g)]< , (5) 
'->l^xp%^ + -"-J 
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where oq = Cp + a^ and Oi = ap — an are combinations of the DM couphngs to protons a.p and 
neutrons a„, and Soo{q), S'oi(g), >S'ii(g) are the spin-dependent nuclear structure functions, 
which we compute according to [39] for ^^^Cs and according to [40] for all other nuclei. In 
the fits we will assume Op = 1, a„ = 0, a choice that leads to a good fit of data. The reason is 
that a spin-dependent interaction couples predominantly to un-paired nucleons. Hence, for 
coupling only to protons DAMA can be compatible with the remaining experiments since 
^^'^I has odd Z, but even N (i.e., one unpaired proton), while ''^Ge, ^^^Xe, ^^^Xe all have 
unpaired neutrons but even Z. For the other extreme choice Op = 0, a„ = 1 the DAMA 
signal is safely incompatible with the other searches, see e.g. [12]. Simultaneous couplings to 
proton and neutrons can still be allowed, depending on the relative size of the two couplings. 

III. DESCRIPTION OF EXPERIMENTS AND ANALYSIS METHODS 

In our fits we use the most sensitive experimental data sets available to date, coming from 
the experiments DAMA/LIBRA, CDMS-II, XENONIO, ZEPLIN-III, CRESST-II, KIMS, 
and PICASSO. In this section we describe each experiment in turn and comment on the 
sensitivity to each of the four classes of DM interactions - eSI, eSD, iSI, iSD. 

A. DAMA 

In 2008, the DAMA collaboration has published results of the combined DAMA/Nal and 
DAMA/LIBRA experiments [1, 2], corresponding to an exposure of 0.82 ton yr for a target 
consisting of radiopure Nal(Tl) crystals. They observe an annual modulation in the signal, 

S{E,t) = So{E) + A{E) cos uj{t- to) , (6) 

where u = 2tt /I yr, to = 152 days. In the fit we use the signal region from 2 to 8 keVee of 
the spectrum from the combined DAMA/Nal and DAMA/LIBRA ("DAMA", for brevity) 
data given in fig. 9 of [2], divided into 12 bins. The data points above 8 keVee are consistent 
with no modulation. Since no signal is predicted in that range from our DM models, they 
do not provide an additional constraint on the fit, and are thus ignored. Our analysis of the 
DAMA data is analogous to the one presented in [11, 29]. 

The signal in DAMA is the energy deposited in scintillation light, while the scattered 
nucleus is loosing energy both electromagnetically and through nuclear interactions (phonon 
excitations). This effect is taken into account by the quenching factors that convert the total 
nuclear recoil energy Ed to the energy seen in the event by the experiment, q x Ed, with units 
of equivalent electron energy (keVee). For DAMA gNa = 0.3 and qi ~ 0.09 [41]. It is, however, 
known [42, 43] that some recoil nuclei, namely those travelling along the crystal planes, 
will not suffer from quenching, but deposit essentially all their energy electromagnetically 
(corresponding to q = 1). The fraction of recoil nuclei for which this happens has been 
calculated in [43], but this calculation leaves room for some debate [14, 44, 45], since the 
channeling effect has not been measured in the relevant energy range so far. In the following 
we will discuss to what extent possible explanations of the global data rely on the presence 
of the channeling effect. 
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For the fit to the DAMA data we construct a function from the annually modulated 
part of the signal 



xLmaK.c^.) = E (^ ^-^"^"^^;^^^ ~ j , (7) 

where the sum is over the energy bins, and A°^'^ (cTj) are the experimental data points (errors) 
in figure 9 of [2] . We also impose a constraint that the predicted unmodulated signal 5*0 from 
DM scattering should not exceed the one measured by DAMA (which consists of background 
and signal) for any given energy bin. A recent simulation provides a background estimate 
for DAMA [46] . While this information could be useful for future more refined fits, we take 
at this point a conservative approach and let the background float freely in the fit. 

We find the best fit point by minimising eq. (7) with respect to WIMP parameters. 
Allowed regions in the (m^,(jp) plane for elastic scattering or in the {m^,ap,S) space for 
inealstic scattering at a given CL are obtained by looking for the contours x^ij^x^'^p) ~ 
Xmin + ^X^(CL), where Ax^(CL) is evaluated for the corresponding degrees of freedom 
(dof), e.g., Ax^(90%) = 4.6 or Ax^(99.73%) = 11.8 for 2 dof. 



B. CDMS-II 

The most recent analysis of CDMS-11 was performed on data taken between July 2007 
and September 2008 in four periods. Only Ge detectors were used for the DM search with a 
total exposure of 612 kg days. Two events were seen in the 10-100 keV energy window, with 
recoil energies of 12.3 keVnr and 15.5 keVnr\ while 0.8 ± 0.1 ± 0.2, 0.04l|J:[]3 and 0.03-0.06 
background events are expected from misidentified surface events, cosmogenic background 
and neutron contamination, respectively. In our fits we take into account that the signal 
efficiency drops from 32% at 20 keVnr to 25% at both 10 keVnr and 100 keVnr by linear 
extrapolation. We also include the previous CDMS search with null result for exposure of 
397.8 kg days obtained between October 2006 and July 2007 [47]. We use a constant energy 
resolution of 0.2 keV for the CDMS germanium detectors. 

In most part of our work we follow the CDMS collaboration and use the data only to set 
an upper bound on a possible signal from DM. To this aim we employ Yellin's maximum gap 
method [48], which by construction leads only to a bound (and never to a positive signal), 
without any assumptions on the possible origin (background or signal) of observed events. 
Only in sec. IV B we are more speculative, and perform a maximum likelihood fit to the two 
observed events assuming a model for the background (details of that analysis are given in 
sec. IV B). 



Here, keVnr refers to the actual nuclear recoil energy, as opposed to the equivalent electron energy reported 
by DAMA. 
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C. XENONIO 



The XENONIO experiment ("XENON", for short) searches for DM scattering on xenon 
nuclei by measuring simuhaneously the scintillation and ionization signals in purified liquid 
xenon. Using a fiducial mass of 5.4 kg they collected a data sample of 316.4 kg day between 
October 6, 2006 and February 14, 2007 [4]. The latest analysis [5] of this data yields 13 
events with an expected background of 7.4 in the 2.0-75.0 keVnr window. This is a re- 
analysis of the same data used in [4], where 10 events had been found in a smaller energy 
window from 4.5-26.9 keVnr. There is still some controversy concerning the effective light 
yield in liquid xenon, which is needed to translate the observed ionization signal into a 
nuclear recoil energy [49, 50]. In our analysis, we use the correction factors from [49], but 
we also discuss the impact of using instead the data from [50]. For the detection efficiency 
and background estimates, we use the numbers given in table I of [5]. The energy resolution 
is computed according to [12] 



with a = 0.579 and /3 = 0.021. We analyze the XENON data using the maximum gap 
method [48]. 

D. ZEPLIN-III 

Like XENONIO, the ZEPLIN-III experiment is also a two-phase liquid xenon time projec- 
tion chamber experiment that has accumulated an exposure of 847 kg days between February 
27th and May 20th 2008 [51]. They observed 7 events in the 2-16 keVee energy window 
(shown in fig. 16 of [51]), which corresponds to the nuclear recoil energy window of 10.7- 
30.2 keVnr. In our analysis we convert keVee into keVnr (and vice-versa) using the scintil- 
lation light yields from fig. 15 of [51]. The detection efficiency is taken from fig. 14 of [51], 
and the energy resolution is assumed to be the same as in XENONIO, eq. (8). We analyze 
the data by employing the maximum gap method. 

E. CRESST-II 

The commissioning run of CRESST-II [19] from March 27th to July 23rd, 2007 has 
a cumulative exposure of 47.9 kg days on CaW04 crystals. Following [19], we use only 
tungsten recoils in our analysis. There are 3 events in the 10-40 keVnr energy window at 
about 17 keVnr, 18 keVnr, and 33 keVnr. We assume a constant efficiency of 90% and a 
constant energy resolution of 1 keV [25]. For the statistical analysis, we again employ the 
maximum gap method. 



The KIMS experiment has an exposure of 3409 kg days taken with low background CsI(Tl) 
crystals [16]. Since no information on the actual running periods was available to us, we 




(8) 



F. 



KIMS 
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assume them to be equally distributed throughout the year. We compute the energy reso- 
lution according to eq. (8) with a = 0.582 and (3 = 0.0021 [52] and assume the detection 
efficiency to be 30% at 3 keVee and 60% at 5 keVee [53]. Between these two energies, we 
use linear interpolation. The quenching factor is computed from a fit to fig. 4 of [52]. The 
KIMS data is reported in a window ranging from 3 keVee to 11 keVee, divided into eight 
bins for each of the four detector modules. We use the Ax^ method to analyze this data in 
the same way as the unmodulated DAMA data. To be conservative, we add penalties to the 
only where the predicted event rate is larger than the measured one, in order to make 
sure that only an upper bound is obtained. 

G. PICASSO 

The PICASSO experiment at SNOLAB [17] is very different from the other experiments 
considered in this work because it is based on the superheated droplet (bubble chamber) 
technique to search for DM recoiling on ^^F nuclei in a C4F10 target. As we will see below, 
scattering on ^^F is very sensitive to SD interactions. The experimental procedure is to 
measure the bubble formation rate as a function of the temperature T. Since bubble forma- 
tion is possible only above a certain threshold energy E'thr which depends on T and on the 
particle species, a DM signal would manifest itself as an increase in the bubble formation 
rate over a certain T interval. We compute E'thr according to eq. (2) of [17]. The PICASSO 
data, corresponding to a total exposure of 13.75 ± 0.48 kg days, is analyzed separately for 
the two best detectors (modules 71 and 72) with the Ax^ method, where the background is 
assumed to have the form a[l + tanh b{T — To)]/2.0, and a, b, and Tq are free parameters 
determined by the fit. We also allow for a 14% uncertainty in the predicted signal [17]. 

IV. RESULTS 

In this section we present the results of our fits to the data. We discuss elastic SI and SD 
scattering in sec. IV A, and inelastic SI and SD scattering in sec. IV C. While in those two 
subsections we follow the CDMS collaboration and use their latest data only to set upper 
limits, in sec. IV B we become more speculative and discuss a possible interpretation of the 
two observed events ClS db SI gnal from SI or SD elastic DM scattering. 

A. Elastic SI and SD scattering 

In fig. 1, we summarize our results for elastic scattering. We show the DAMA allowed 
regions compared to constraints from other experiments in the plane of DM mass and 
the interaction cross section for SI scattering and SD scattering off protons. We observe that 
in both cases the DAMA regions are excluded by the bounds. For the SI case (left panel) 
the most important bounds come from CDMS and especially from XENON. They exclude 
the DAMA regions regardless of the assumptions on channeling at high CL. Let us note that 
for the DAMA region around 10 GeV without channeling additional contraints from CDMS 
data on silicon apply (not shown) [11]. 
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FIG. 1: DAMA allowed regions (90% and 3(j CL) and constraints from other experiments 
(90% CL) for SI scattering (left) and SD scattering off protons (right). Shaded DAMA regions 
have been obtained assuming the channeling effect according to [43] , while the black contour curves 
correspond to no channeling. 

To study the impact of the new 2009 analyses from XENON [5] and CDMS [3], we compare 
in fig. 2 the old and new data sets, and show also the impact of different assumptions on 
the effective light yield Lgfj in XENON. We observe that despite the large exposure, the 
new CDMS data has very small impact on the lower bound on m^. The reason is that the 
two observed events are located at low energies, which are most relevant if is small, and 
therefore the maximum gap method leads to a not so strong limit. The new CDMS data is 
more important if is large so that a signal is expected also at larger recoil energies, where 
no events have been observed. Therefore, thanks to the large exposure the limit improves. 
We also note a rather significant improvement in the low-mass limit from XENON due to 
the 2009 analysis. There are two reasons for this effect: first, the energy threshold has 
been lowered from 4.6 keVnr [4] to 2 keVnr [5], and second, the one event located close to 
4.6 keVnr in the 2007 analysis has been eliminated in the 2009 analysis. This leads to a 
large energy interval at low recoil energies without events, which improves the limit for low 
masses. 

Considering the case of SD elastic DM-nucleus scattering (right panel of fig. 1), we 
observe that rather strong constraints come from PICASSO, if DM couples mainly to protons 
(ttp = 1, a„ = 0). Assuming channeling according to [43] the DAMA region at 90% CL is 
excluded by the 90% CL bound from PICASSO, while both experiments are marginally 
compatible at 3(j. Without the channeling effect, there is no overlap of allowed regions. We 
do not show the case of SD scattering off neutrons (op = 0, a,i = 1), since in that case the 
DAMA region is safely excluded by CDMS and XENON [12], see also fig. 3. The reason is 
that the ^^F nuclei in the PICASSO experiment have an unpaired proton, while the spin- 
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FIG. 2: Comparison of different data sets for CDMS and XENON for SI elastic scattering. 
XENON 07 refers to the analysis from [4], while XENON 09 corresponds to the re-analysis of 
the same data from [5]. The CDMS 2008 and 2009 data sets are from [47] and [3], respectively. For 
XENON, different assumptions on the effective light yield Leg are used. Thick curves are based 
on the measurement [49], while thin curves are based on the alternative data set [50]. 

sensitive isotopes in CDMS (^^Ge) and XENON (^^^Xe, ^^^Xe) have an unpaired neutron. 
Let us mention that in the proton case, also the COUPP [15] experiment provides a relevant 
constraint. However, it was not possible for us to implement a simulation of COUPP using 
the available information. 



B. A signal in CDMS? 

Our default analysis of CDMS data uses the maximum gap method [48], which produces 
by construction only an upper limit on a DM signal. Let us now be somewhat more specu- 
lative and interpret the two events observed in the latest CDMS data as positive signal from 
DM scattering (see also [54]). In order to do this we use a model for the energy shape of the 
expected background from surface events. According to [3], this background is estimated by 
using the measured events in the 2o window from their previous analysis, where no timing 
cut is yet imposed, see figure 3 of [47]. We count the number of events in the signal region 
and perform a fit to the energy distribution. Normalizing to the expected total number of 
background events (i.e. 0.8), we find dN"^^^' / dE^ = -0.00295 + 0.463/^^ where the recoil 
energy is in keVnr. Using this parameterization for the expected background shape we 
perform a fit to the two events by using the so-called extended maximum likelihood method 
[55]. 

In fig. 3 we show the CDMS allowed regions compared to the constraints from the other 
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FIG. 3: Allowed regions for CDMS 2009 data {la, 90% and 3a CL), DAMA (90% and 3a CL), 
and constraints from other experiments (90% CL) for elastic SI scattering (left), SD scattering off 
protons (middle), and SD scattering off neutrons (right). 

experiments for SI and SD elastic scattering. CDMS allowed regions are defined by Ax^ = 
— 21og£/£™^^ contours for 2 dof. We find that at la CL a closed allowed region appears 
for CDMS ("positive signal"), while already at 90% CL only an upper bound is obtained. 
The CDMS favoured region is largely excluded by the XENON bound (and the PICASSO 
bound, in case of SD scattering off protons). In the case of SI scattering we have checked 
that in the combined fit to all all experiments except DAMA no closed region appears even 
at 1(7, and we obtain only an upper limit. 

We have also verified that in the case of inelastic scattering no closed region appears for 
CDMS. This follows from the fact that for inelastic scattering the signal is shifted to larger 
recoil energies, where CDMS sees no events, and therefore the data do not favour a positive 
signal. Hence, in the following we return to the conservative approach and use CDMS data 
only to set an upper limit on a DM signal using the maximum gap method. 



C. Inelastic SI and SD scattering 

Let us now consider the assumption of inelastic DM scattering, put forward in ref. [18] in 
order to reconsile the DAMA signal with constraints from other experiments exploring the 
modified kinematics, see eq. (2) and the discussion given there. We start by first presenting 
results for the SI case, largely discussed in the literature, and then extend the inelastic 
scattering hypothesis also to SD interactions. 

In order to identify possible solutions we have performed a combined analysis of all ex- 
periments, by adding the functions of the individual experiments. For those experiments 
which we analyse with the help of the maximum gap method (see sec. Ill) we proceed as 
follows. For a given point in the parameter space the probability obtained by the maxi- 
mum gap method is converted into a Ax^ by inverting the integral over the x^-distribution 
for 2 dof. This "fake" Ax^ is added to the one from the remaining experiments. Then 
we perform a scan over all three model parameters, m^, ap, and 6, in order to search for 
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local minima in the total x^- While this recipe to incorporate the maximum gap method 
leads only to approximate confidence regions, it suffices for our purpose to locate potential 
solutions in the 3-dimensional parameter space. 

Fig. 4 shows the results of such a scan for the SI case. We show projections of the 3- 
dimensional allowed regions onto the three 2-dimensional planes. We identify three local 
minima. The global minimum is Xmin.giob — ^O-O located at 

(Tp = 2 X 10"^^ cm^ , m^ = 10.1GeV, (5 = 39keV. (9) 

This solution corresponds to channeled events from iodine [25] and relies on the channeling 
calculations from [43] (scattering on sodium does not contribute to the signal, since the 
kinematics of inelastic scattering favor heavy over light nuclei). Furthermore this solution 
corresponds actually to a tiny, rather fine tuned region in 6 and m^. The values are choosen 
such that the minimal DM velocity fmm required to give a recoil energy of 3 keV needed to 
explain the DAMA signal is very close the galactic escape velocity. Indeed, the parameters 
are tuned such that the signal in DAMA is non-zero only in summer but zero in winter. 
This maximally enhances the modulated signal, while at the same time suppresses the 
unmodulated rate as well as the signal in the other experiments.^ This is illustrated in the 
upper panel of fig. 4, where we show that the allowed region is located precisely between the 
contour curves for Vmin = Vesc in summer and in winter. Because of this fine tuning the best 
fit value for the cross section given in eq. (9) is very sensitive to the precise implementation of 
the DM halo profile and minor modifications in the analysis. This tuning of DM parameters 
S and relative to properties of the galactic halo (fesc) is rather un-natural, and we 
consider this solution as being disfavoured despite the formally very good value. The same 
arguments apply for a local minimum around ap = 10^'^^cm^,m^ = 40GeV,5 = 130 keV, 
corresponding to quenched events on iodine. 

Another local minimum appears around ap = 10~^^cm^,m^ = 50 GeV, 5 = 130 keV 
(light shaded regions in fig. 4). This corresponds to the conventional iDM solution discussed 
recently by many authors, e.g. [21-26]. This solution does not suffer from the ffne tuning 
problem, but we find that it is disfavored with respect to the best fit with Ax^ = 16.1. The 
reason are strong constraints mainly from CRESST, which is optimal for inelastic scattering 
since the heavy tungsten target nuclei allow for very efficient DM energy loss. In order to 
illustrate these constraints we show in fig. 5 allowed regions in the ap,m^ plane for fixed 
values of S. The upper and lower left panels correspond to the SI fits with S = 130 keV and 
40 keV, respectively. We observe that the 3a DAMA region at 5 = 130 keV is completely 
excluded by the CRESST 90% CL bound and quite strongly constrained by several other 
experiments. The very thin strip corresponding to the DAMA region in the lower left plot 
illustrates the fine tuning problem. Furthermore, we observe that this solution is marginally 
compatible with the XENON bound. 

Let us now move to the spin-dependent case. The results of the parameter scan of the 
combined is shown in fig. 6. We observe that the allowed regions are much larger than 



^ We remark that in such a case the fit to the DAMA modulated signal in terms of a cosine according to 
eq. (6) might not provide a good description, since the signal would correspond to a truncated cosine. We 
are not exploring such additional signatures here. 
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FIG. 4: Global fit of inelastic, spin-independent DM. We show projections of the 3-dimensional 
regions at 90% and 3ct CL onto the three 2-dimensional planes, by minimizing the global in 
each case with respect to the third (un-displayed) parameter. Confidence regions are defined for 
2 dof. The fit includes CDMS (2008 + 2009 data), XENON (2009 analysis), DAMA, CRESST-II, 
ZEPLIN-III, and KIMS. The dark shaded regions are the allowed regions defined with respect 
to the global minimum, while the light shaded regions are defined relative to a local minimum, 
which by itself is disfavored relative to the global minimum with Ax^ = 16.1. The open contours 
correspond to DAMA data only. In the upper panel we show also contours of Vyam = ^^esc + 'dearth 
for a recoil energy oi E = 2> keV, with and without quenching of iodine scatters, where dearth is the 
velocity of the earth relative to the halo, depending on the time in the year. 
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FIG. 5: DAMA allowed regions (90% and 3a CL) and constraints from other experiments 
(90% CL) for inelastic DM scattering with SI interactions (left) and SD interactions with pro- 
tons (right). We show regions in the {m^,ap) plane (2 dof) for fixed DM mass splitting 6. The 
upper panels {5 = 130 keV) correspond to a signal in DAMA from quenched events on iodine, 
whereas the lower panels {6 = 40 keV) correspond to channeled events on iodine according to [43]. 

in the SI case. There are two different regions corresponding to fits of similar quality, 
corresponding to quenched events in DAMA and channeled events according to [43]: 

~ 40 - 70 GeV, 6 ~ 130 keV (quenched events) (10) 
~ 10 GeV, S ~ 40keV (channeled events) (11) 
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niy [GeV] 
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FIG. 6: Global fit of inelastic, spin-dependent DM. We show projections of the 3-dimensional 
regions at 90% and 3ct CL onto the three 2-dimensional planes, by minimizing the global in 
each case with respect to the third (mi-displayed) parameter. Confidence regions are defined for 
2 dof. The fit includes CDMS (2008 + 2009 data), XENON (2009 analysis), DAMA, CRESST- 
II, ZEPLIN-III, and KIMS. (We omit PICASSO which, due to the light target nucleus ^^F has 
virtually no sensitivity to inelastic DM). The shaded regions refer to the global fit, whereas the 
open contours correspond to DAMA data only. In the upper panel we show also contours of 
Vmm = 'Wesc + dearth for a recoil energy of -E = 3 keV, with and without quenching of iodine scatters, 
where dearth is the velocity of the earth relative to the halo, depending on the time in the year. 
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FIG. 7: projected onto the ap,m^, and 6 axes, where we minimize with respect to the two 
un-displayed parameters. We show the for DAMA data only (orange), and for the global data 
(black) assuming iSI interactions (solid) and iSD interactions on protons (dashed). 

with cross sections in a wide range of few x 10""^^ cm^ ^ CTp ^ few x 10~^^ cm^. The right 
panels of fig. 5 show that the allowed regions are safely compatible with the constraints 
from all other experiments. The explanation is as follows: the SD coupling to protons dras- 
tically reduces the power of even Z target experiments (XENON, CDMS), while the inelastic 
kinematics strongly disfavour light targets (PICASSO), which provide a main challenge for 
DAMA in the elastic SD case, see fig. 1. 

It is, however, important to remark that for each fixed 5, it is only a very small range 
of that gives a good fit to the DAMA data. Varying by a small amount requires 
a large change of the cross section to maintain a good fit, see fig. 5 right panels. This is 
because the scattering is sensitive to the exponentially suppressed tail of the DM velocity 
distribution close to fesc- However, in contrast to the SI case, we observe from the top panel 
of fig. 6 the 90% CL regions extend relatively far way from the v^nm = v^sc curves. Hence, 
the iSD case does require a relatively precise tuning of model parameters {ap,m^,6) among 
themselves, but the tuning with respect to astrophysics is not necessary here. Therefore, if 
we assume the spin-dependent inelastic scenario to be true, the results can be interpreted 
as reflecting very high sensitivity of direct detection experiments to the model parameters 
up to astrophysical uncertainties. 

Let us mention that we have not been able to include data from the COUPP [15] ex- 
periment in our analysis, due to missing information. In contrast to the elastic SD case we 
do not expect that COUPP will provide a relevant constraint in the inelastic case. Namely, 
the scattering on ^^F is negligibe due to its small mass, in the same way as for PICASSO. 
The constraint from iodine contained in their CF3I target should also be much weaker than 
the one coming from KIMS due to the much larger exposure of the latter (3409 kg days 
for KIMS vs. 250 kg days for COUPP). We do not show the results for iSD scattering off 
neutrons, but we have checked that for the same reasons as in the elastic case, the DAMA 
region is safely excluded by constraints from XENON and CDMS. 

To conclude this section, we show in fig. 7 the projections of the DAMA-only and global 

function for the three parameters dp, m^, 6, separately, for the iSI as well as iSD scenarios. 
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V. A SIMPLE MODEL 



As we have seen in the previous section, the possibihty that DM scatters inelastically 
with spin-dependent cross section provides a viable explanation for the DAMA signal that 
simultaneously avoids the bounds from the other direct detection experiments. We now 
show that iSD scattering can be realized in concrete models. 

It actually does not require too much work to find such a model. To be concrete let 
us start with fermionic DM that interacts with the visible matter through the effective 
interaction ^ 

-^int = [^roMV^] [q^visq] , (12) 

where Fdm.vis denote the Dirac structure of the four-fermion operator which we keep general 
for now, ip = (?7,C^) is a Dirac fermion (for two-component spinors we use the notation 
of [56]) and q are the light quark fields. We could have equally well chosen couplings to 
leptons, in which case the analysis of DAMA would have followed ref. [29]. The above 
effective interaction can arise from an exchange of heavy mediators with mass 0{A) under 
which both visible and DM fermions are charged. In [18] the four-fermion interaction (12) 
was chosen to be of the V<^V form (i.e., Fdm = Tvis = l'^)- This gives a spin-independent 
scattering cross section for inelastic DM. We show below that tensor interactions give spin- 
dependent inelastic scattering instead. To the best of our knowledge this realization has not 
been discussed in the context of inelastic DM scattering in the literature before (for some 
particle physics realizations of inelastic DM see [23, 24, 57-62]). 

If ijj is a. Dirac fermion, the interaction (12) leads to elastic scattering. However, if in 
addition to the Dirac mass term miljilj there are also Majorana mass terms {Srjrjri + (5g^^)/2, 
then the Dirac fermion splits into two Majorana fermions with masses m±6 (for simplicity 
let us take 6^^ = 6^ = 6). The mass eigenstates are [18] 

Xi = ^{v-0, ^2 = ^(r? + 0- (13) 

It is reasonable to assume that 6 m, since Majorana mass terms break a global symmetry, 
while the Dirac mass term conserves it. This is exactly the hierarchy needed phenomeno- 
logically, since m ~ C(100 GeV), 6 ~ C(100 keV) are needed for DAMA. 

Let us now assume that the four-fermion interaction between DM and the visible sector 
(12) is of T ® T form, 

£int = ^ [^E,,V^] [qE^^q] , (14) 

where T,^'^ = From the relations Xi'^^'^Xj = ~Xj(^^'^Xi: xl(^^'^x] = ^Xj^^'^x! 

(see e.g. [56]), it follows that for Majorana fermions the diagonal tensor operator vanishes.^ 
Thus, one finds for the DM tensor current 



3 Here, a^"' = iia^'a" - a''a'')/4:, af" = {{af^cr" - 5-''cr^)/4, and ct^ = {l,cr'), ct^ = (l,-o-'), with cr* the 
Pauli matrices. 
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which leads to inelastic scattering for 5 7^ 0. Furthermore, it is well known that in the 
nonrelativistic limit the T ®T interaction leads to spin dependent scattering, see e.g. [63]. 
For instance, in our case the matrix element for Xi ~^ X2 scattering on a nucleus N has 
the form {N\qSiq\N){(l ^,ai(i^s), with Si the spin operator and C2,s',Ci,s the nonrelativistic 
two-component spinor wave functions for the DM particles. 

It is important to note that the T ®T current is the only chiral structure that leads to 
spin-dependent inelastic scattering in the above simple model of two Majorana fermions split 
by small Majorana mass terms. The T ®TA structure of the four-fermion interaction would 
vanish in the f — )■ limit, while the TA ® TA interaction is equivalent to T ® T as is easily 
checked from the definition of 75. The V ®V product leads to SI inelastic scattering [18], 
A® A to SD elastic scattering, the V ®A product vanishes in the nonrelativistic limit, while 
the scalar and pseudoscalar couplings obviously do not lead to spin dependent interactions. 

One still has the freedom to choose appropriate values for Wilson coefficients Ct for 
different flavors of the quark current, i.e. the couplings to u and d quarks. In section IV we 
explored two extreme cases where the coupling is proportional to the charge (so that DM 
scatters only on protons) or that the coupling to d is twice as large and opposite to the 
coupling to u quarks (so that DM scatters only on neutrons). From the flts the flrst option 
is preferred. 

One can also construct a model with DM that is a scalar and that has inelastic spin- 
dependent scattering on visible matter. But this interaction is inevitably suppressed by 
the nonrelativistic DM velocity because of the derivative in the DM current and thus less 
realistic phenomenologically in light of relatively large scattering cross sections needed to flt 
DAMA. 

VI. CONCLUSIONS 

The CDMS-II collaboration has recently reported the observation of two events, with 
an expected background of 0.9 events [3], by about doubling the exposure with respect to 
previous results. We have performed a combined analysis of direct DM searches including 
the claimed signal by DAMA [1, 2], the recent CDMS-II results, as well as a re-analysis 
of XENONIO data. We have considered four classes of possible WIMP-nucleus scattering 
models: elastic (e) or inelastic (i) scattering and spin-dependent (SD) or spin-independent 
(SI) scattering. This covers a large set of DM models. While three types of WIMP-nucleus 
scattering have already been considered in the literature (i.e., eSD, eSI, and iSI scattering), 
our analysis is the flrst to also include the fourth possibility — the inelastic spin dependent 
scattering. In fact it is this latter possibility that can simultaneously explain the DAMA 
signal and avoid bounds from the other direct detection experiments as we demonstrated in 
the present paper. 

For eSI scattering the DAMA region is safely excluded by data from XENONIO, based 
on a recent re-analysis of their data [5], mainly due to a lower energy threshold, while for 
the eSD case DAMA is in conflict with the bound from PICASSO. For iSI scattering we flnd 
two possible solutions. The traditional "iDM" region with ~ 50 GeV and 5 ~ 130 keV 
is disfavoured by CRESST-II data and further constrained by CDMS-II. A low mass region 
exists around ~ 10 GeV and (5 ~ 40 keV. It requires, however, the presence of the 
channeling effect in DAMA, and leads to rather severe tuning of the DM mass and mass- 
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splitting with respect to the properties of the galactic DM halo, such that the WIMP velocity 
distribution is sampled precisely around the escape velocity. 

For the iSD case, on the other hand, we do find allowed regions of the parameter space 
such that the signal in DAMA is explained, while there is no confiict with any of the other 
experiments for non-vanishing couplings to protons and suppressed couplings to neutrons. 
The solution suffers from some tuning among the DM parameters ap,m^, and 6 (which 
a priori is not a problem, but just indicates very high sensitivity of the experiments to 
these parameters), while they do not require a very precise tuning with respect to the 
escape velocity. The required DM parameters for the two solutions are around 40- 
70 GeV, 6 ^ 130 keV (corresponding to quenched events in DAMA), and ~ 10 GeV, 
5 ~ 40 keV (which relies on the channeling effect). Cross sections are found in a wide range 
of few X 10~^^ cm^ ^ < few x 10~^^ cm^. We presented also a simple toy model, showing 
how to realize iSD scattering in a specific framework. 

In the case of elastic scattering, we also adopted a somewhat speculative approach to 
the recent CDMS-II data, and have performed a maximum likelihood fit to the two events. 
We find that at Icr CL a closed allowed region appears for CDMS ("positive signal"), while 
already at 90% CL only an upper bound is obtained. The CDMS favoured region is largely 
excluded by the XENON bound (and the PICASSO bound, in case of SD scattering off 
protons). In a combined analysis of CDMS and the exclusion limits from other experiments, 
the "signal" becomes less than la. More information on this can be expected soon from 
XENON- 100 [64]. 

Note added in proofs 

After this work was completed, we were able to perform a fit to the data from the COUPP 
experiment. This was made possible by kind assistance from the COUPP collaboration, 
especially Juan I. Collar, who provided crucial data on the bubble formation threshold 
energies. We have checked that COUPP limits are subdominant in the case of eSI, iSI, and 
iSD scattering. For eSD scattering, COUPP improves the PICASSO limit by up to a factor 
of 2 in (Tp for > 50 GeV. 
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FIG. 8: CoGeNT and DAMA allowed regions (90% and 3cr CL) and constraints from other ex- 
periments (90% CL) for elastic spin-independent scattering. Shaded DAMA regions have been 
obtained assmning channeling according to [43], while the black contours correspond to no chan- 
neling. The CoGeNT upper bound follows from a fit including an exponential background, while 
the allowed region is obtained when the exponential is removed from the background model. 

Appendix A: Interpretation of CoGeNT results 

After publication of this manuscript, the CoGeNT collaboration has published results 
from the initial run of their ultra low noise germanium detector in the Soudan Underground 
Laboratory [65]. The detector has a very low energy threshold of 0.4 keVee, the lowest 
achieved so far by any dark matter experiment. The data from an eight week exposure 
of the detector with a fiducial mass of 330 g reveals ~ 100 events near the low-energy 
threshold for which the collaboration was not able to identify a background source. The 
"signal" is consistent with an exponential background. It is, however, also compatible with 
the hypothesis of an (9(10 GeV) WIMP. Below we show that this hypothesis is strongly 
disfavored by other experiments, though not completely ruled out. 

In the analysis we assume standard astrophysical parameters for the dark matter halo, 
see Section II. In the fit we use the data from Fig. 3 of [65] in the 0.4-3.2 keVee energy 
range. For the quenching factors we use the approximate formula (based on [66], Eq. (2.6), 
with K = 0.2 [67]) 

, ^ , (Al) 

1 + ^l + lh.hh/E 

where E is the energy in keVee. Whether or not one finds in the fit a DM signal depends cru- 
cially on the assumed background. We perform fits for two different background models: (i) 
a 7-parameter model consisting of a constant (1 parameter), an exponential (2 parameters), 
and two Gaussians at the known positions of the ^^Ge and ^^Zn lines (2x2 parameters); 
this background model is similar to the one used by the CoGeNT collaboration, (ii) A 
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FIG. 9: CoGeNT and DAMA allowed regions (90% and 3a CL) and constraints from other 
experiments (90% CL) for SD scattering on protons (left) and on neutrons (right). See also caption 
of Fig. 8. 

similar background model, but without the exponential. The results are shown in Fig. 8 for 
elastic SI scattering, and in Fig. 9 for elastic SD scattering. We also include fits to the data 
from the CDMS Sihcon detectors (labeled CDMS-Si) [68, 69]. In the SD case, we use the 
structure functions for ^^Si as given in [40]. 

For the full 7-parameter background model only an upper limit on the WIMP cross-section 
is obtained for both SI and SD scattering. It is only if one assumes that the background 
does not contain an exponential component that a closed allowed region appears in the 
— ap plane, with ~ 5 — 10 GeV, cTp ~ 10~^ pb for SI scattering, and ~ lO'' pb 
(cr„ ~ 10 pb) for SD scattering on protons (neutrons). In the case of SD scattering on 
neutrons, the CoGeNT-preferred region overlaps with the DAMA-preferred region. For SI 
scattering, marginal overlap between CoGeNT and DAMA could be achieved if the frac- 
tion of channeled events is assumed to be smaller than what was assumed by the DAMA 
collaboration (orange shaded contours), but larger than zero (black contours). Also, an 
admixture of an exponential background to a DM signal in CoGeNT could shift the allowed 
region to lower values of a^, closer to the DAMA allowed region. However, without a clear 
prediction for the background all such modifications are mere speculations. Furthermore, 
both the DAMA allowed region and the CoGeNT "signal" region are ruled out at 90% CL 
by the other experiments for the standard choices of experimental parameters. XENON-10 
and CDMS-Si rule out the SI case as well as SD scattering on neutrons, while COUPP and 
PICASSO rule out SD scattering on protons as an explanation for CoGeNT and/or DAMA, 
cf. Figs. 8 and 9. As fig. 2 shows, relaxing the assumptions on the effective light yield L^s 
in XENON-10 can make XENON-10 and CoGeNT (but not DAMA) marginally compatible 
at 90% C.L. We have also checked that CoGeNT results do not further constrain inelastic 
DM (whether SI or SD) discussed in Section IV C. 

In conclusion, we find that the background only hypothesis gives an excellent fit to the 
CoGeNT data, while a DM interpretation of this data is disfavored by other experiments. 
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We differ in this conclusion from the authors of [70], who find that CoGeNT and DAMA 
are consistent with the remaining experiments. This difference may be traced back to 
the treatment of L^g for XENON-10. The authors of [70] assume the true L^g to he at 
the lower end of the 1-a error bars of the most conservative measurement available [50]. 
Given the discrepancy between different measurements of L^s, the possibility of such a 
systematic shift cannot be excluded until systematical errors in the L^s measurement are 
better understood. Other differences between our analysis and that of [70] are in the choice 
of DM halo parameters and in that the authors of [70] demand that the background does 
not exceed the DM signal, while we allow both to float freely. 
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